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INVESTIGATION  OF  THE  ACOUSTIC 
EMISSION  AND  DEFORMATION 
RESPONSE  OF  FINITE  ICE  PLATES 

P.C.  Xirouchakis,  M.  Chaplin  and  W.F.  St.  Lawrence 


INTRODUCTION 

This  report  is  the  result  ol‘ a  preliminary  investiga¬ 
tion  undertaken  to  examine  the  mechanical  properties 
of  ice  undergoing  flexural  deformations.  The  primary 
purpose  of  tilts  work  was  to  study  the  time-dependent 
deformation  of  small  simply  supported  ice  plates  and 
the  microfracture  properties  as  indicated  by  acoustic 
emission  activity  .  Since  the  experiments  conducted 
generally  led  to  catastrophic  failure  of  the  ice  plates, 
we  were  able  to  obtain  information  both  on  the  loads 
that  caused  the  plates  to  fracture  and  on  the  time  to 
failure  of  the  ice  plate.  Although  this  work  represents 
a  preliminary  effort,  we  (eel  it  is  of  some  benefit  in 
terms  of  establishing  a  method  of  testing  finite  ice 
plates.  It  also  gives  us  some  indication  of  magnitude 
of  the  loads  which  induce  significant  time-dependent 
deformation,  microfracture  and  catastrophic  failure 
of  the  ice. 

A  major  problem  encountered  in  our  work  was  re¬ 
lated  to  the  development  of  a  methodology  for  test¬ 
ing  small  ice  plates  subjected  to  loads  that  induced 
flexural  deformation.  Although  a  significant  amount 
of  work  has  been  carried  out  on  columnar-grained  ice 
subjected  to  uniaxial  compressive  loads  (for  example, 
see  Ciold  1973  and  Sinha  I97K)  we  were  unable  to 
find  work  in  the  literature  which  related  to  flexural 
loading  of  small  ice  plates.  Erederking  and  Gold 
( 1974)  did  work  on  edge  loaded  ice  plates  and  the  re¬ 
port  of  their  methodology  was  helpful  in  conducting 
this  work. 

Throughout  the  course  of  our  three-week  investiga¬ 
tion  we  continually  updated  our  experimental  tech¬ 
nique.  We  now  feel  that  a  viable  method  for  examin¬ 


ing  the  flexural  response  of  simply  supported  ice  plates 
has  been  established.  Since  our  specimen  size  was 
small  (plate  surface  area  on  the  order  of  lm2 )  we  were 
not  able  to  load  the  ice  plates  while  they  were  floating 
in  the  water  in  which  they  grew.  Eventually  we  devel¬ 
oped  a  support  system  that  allowed  us  to  measure 
considerable  deflections  of  the  ice  plates  while  also 
allowing  us  to  assume  the  plate  to  be  simply  supported. 

The  elastic  solution  for  the  deformation  of  a  simply 
supported  ice  plate  with  the  configuration  of  the  ones 
tested  is  presented  in  Appendix  A  of  this  report.  By 
applying  the  correspondence  principle  of  linear  visco¬ 
elasticity  to  the  elastic  solution  for  the  plate  deforma¬ 
tion.  we  have  also  obtained  the  time-dependent  visco¬ 
elastic  solution  for  the  deformation  of  a  simply  sup¬ 
ported  plate  corresponding  to  our  test  configuration. 

In  conjunction  with  this  later  solution  we  have  meas¬ 
ured  the  viscoelastic  parameters,  assuming  the  ice  to 
be  a  Maxwell  fluid. 

A  primary  reason  for  conducting  these  experiments 
was  to  examine  the  acoustic  emissions  from  ice  sub¬ 
jected  to  flexural  loads.  The  microfracture  properties 
of  ice  in  uniaxial  compression  have  been  examined  by 
Gold  ( 1972).  However,  no  work  has  been  reported 
which  describes  the  microfracture  of  ice  in  flexure. 

To  this  end  we  monitored  the  ice  plates  using  high- 
frequency  ( 1 00-  to  300-kHz)  acoustic  emission  mon¬ 
itoring  equipment.  The  results  of  this  monitoring 
showed  that  the  ice  plate  does  emit  slress  waves  that 
can  be  defected  as  acoustic  emissions.  However,  al 
this  juncture  the  results  arc  far  from  being  definitive. 
We  do  feel  that  this  is  a  reasonable  line  of  research  to 
pursue. 

In  Ihe  course  of  our  experiments  we  loaded  all  the 


plates  tested  until  they  tailed  catastrophically .  In 
seveial  ul  the  tests  the  fracture  look  place  immediately 
upon  application  of  the  load.  In  several  other  tests  sig¬ 
nificant  time-dependent  deformation  took  place  before 
li  icture  was  obseived.  I  rom  this  we  have  obtained 
so  ue  mini  illation  on  the  loads  required  io  fracture  a 
su.iply  supported  ice  plate. 

As  a  means  of  considering  the  structure  of  the  ice, 
we  made  thin  sections  of  all  the  ice  plates  tested  after 
failing  took  place.  This  analysis  allows  us  to  consider 
out  load  deformation  response  and  also  the  acoustic 
emission  and  fracture  iespon.se  in  light  of  the  grain 
structure  of  the  material.  We  hope  to  use  these  data 
on  ciystal  structures  in  the  future  not  only  to  correlate 
structure  with  imposed  parameters  but  also  in  control¬ 
ling  the  nature  of  the  ice  structure  being  tested.  In 
this  series  of  experiments  no  attempt  was  made  to  in¬ 
fluence  the  crystal  growth:  thus  the  ice  tested  grew 
under  the  influence  of  the  variations  m  the  coldrooin 
environment  and  as  a  result  the  crystal  structure 
varied  significantly. 

It  was  our  plan  to  analyze  and  compare  the  behav¬ 
ior  of  freshwater  ice.  sodium  chloride  ice  and  urea  ice. 
However,  we  found  that  given  the  time  constraint  im¬ 
posed  on  our  experiments  (3  weeks)  this  was  not  pos¬ 
sible  to  accomplish.  In  the  future  we  hope  to  be  able 
to  complete  this  task.  The  purpose  of  this  report, 
which  is  somewhat  altered  in  terms  of  our  original 
goals,  is  1 )  to  document  our  experimental  effort,  2)  to 
establish  the  analytic  basis  for  examining  plate  defor¬ 
mation.  and  3)  to  obtain  preliminary  information  on 
the  acoustic  emission  and  failure  properties  of  ice 
subjected  to  flexural  loads.  We  hope  that  this  report 
will  be  of  some  value  to  those  who  will  continue  to 
work  along  these  lines. 

EXPERIMENTAL  PROCEDURE 
AND  CONSIDERATIONS 

A  significant  portion  of  the  time  allotted  for  this 
research  endeavor  was  consumed  in  developing  our 
experimental  methodology.  The  instrumentation  for 
measuring  the  deformation  and  acoustic  emissions  of 
the  loaded  ice  plate  followed  standard  procedures. 

The  manner  in  which  the  ice  plate  was  grown  and  the 
plate  support  system  used  required  a  significant  num¬ 
ber  of  trial  and  error  attempts  before  we  were  satisfied 
that  we  were  obtaining  the  desired  test  configuration. 

In  this  section  we  give  a  detailed  description  of  the 
procedures  used  in  growing  and  testing  our  ice  plates, 
and  also  we  describe  the  acoustic  emission  and  load 
deformation  system. 


Growth  of  the  ice  plate 

It  was  our  goal  to  use  relatively  thin  ice  plates 
grown  in  a  stress-free  state.  To  facilitate  this  goal,  a 
square  plywood  box  with  internal  dimensions  of  1.02 
x  1 .02  m  and  0.6  m  deep  was  used.  The  inside  of  the 
box  was  coated  with  an  epoxy  resin  to  make  it  water¬ 
tight.  Around  the  outside  perimeter  of  the  box.  at 
approximately  0.45  m  from  the  bottom,  a  heat  tape 
was  secured  to  the  box.  The  heat  tape  was  in  turn  con¬ 
nected  to  a  voltage  control  device  so  tlial  the  heat  gen¬ 
erated  could  he  regulated.  The  box.  which  was  inside 
a  coldrooin.  was  then  filled  with  distilled  water  to  a 
depth  of  about  0.45  m.  so  that  the  walei  surface  was 
at  the  same  level  as  the  heal  tape.  The  coldrooin  for 
these  tests  was  kept  at  a  nominal  temperature  ol  -I0CC . 

By  trial  and  error. a  power  setting  was  found  for 
the  heat  tape  at  which  the  water  at  the  edge  of  the 
box  would  not  frce/c  while  all  othei  surface  water 
would.  To  ensure  that  only  the  surface  water  of  the 
box  would  Iree/e.  the  exterior  sides  of  the  box  were 
insulated  with  6  in.  of  foil-backed  fiberglass  insulation 
willi  a  thermal  resistance  value  lated  at  R  Id.  We 
found  that  a  thin  sheet  of  ice  would  form  oil  the  sur¬ 
face  of  the  water  when  the  temperature  reached  approxi¬ 
mately  -6°('  just  above  the  water  surface  in  the  -lO^C 
coldrooin.  Alter  a  period  of  some  24  to  36  hours  the 
thickness  of  the  ice  plate  would  be  between  I  and  3 
cm  thick,  at  which  point  we  considered  the  plate  ready 
for  testing.  The  ice  plate  temperatures  were  - 1 0cC 
at  the  time  of  each  lesi. 

In  all  the  tests  conducted  here  we  used  distilled 
water  as  our  fluid  and  did  not  make  any  attempt  to 
seed  the  water  surface  to  obtain  a  uniform  grain  size 
in  our  plates.  However,  it  is  possible  that  nucleation 
of  the  ice  crystals  making  up  the  ice  plate  occurred  dur¬ 
ing  the  defrost  cycle  ol  the  coldrooin.  when  a  consider¬ 
able  numhei  of  nucleation  particles  arc  available  in  the 
an  to  seed  the  ice  plate.  To  establish  more  control 
over  our  crystal  structure  in  the  future,  we  will  use  the 
powder  snow  seeding  technique  described  by  I  reder- 
king  and  Gold  ( I ‘>74 ). 

Throughout  file  growth  period  of  the  plate,  both  the 
water  and  aii  temperature  were  monitored  with  copper- 
constantan  theimocouplcs.  When  the  plate  was  fully 
formed,  a  thermocouple  was  also  frozen  into  the  ice. 

The  ice  plate  of  the  size  grown  was  larger  than  the 
plate  we  desired  to  lest.  To  cut  the  plate  to  the  proper 
dimensions  immersion  heaters  (head  bolt  heaters  for 
automobiles)  were  used.  The  temperature  of  the  heat¬ 
ers  was  controlled  with  a  voltage  control  device  (Variac) 
to  limit  the  heat  ouput.  If  we  failed  to  limit  the  power 
output  of  the  heaters  used  for  ice  cutting,  the  heating 
element  would  invariably  burn  out.  During  the  cutting 
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Figure  1.  Ice  plate  support  structure  with  a  plate  positioned  on  the 
structure.  The  area  of  load  application  is  indicated. 


process  we  noticed  no  indication  of  local  ice  fracture 
near  the  cutting  instrument. 

Support  of  the  ice  plate 

A  considerable  effort  was  made  to  obtain  a  simple 
support  system  on  which  the  ice  plate  could  he  tested. 
After  some  trial  and  erroi  we  arrived  at  a  design  which 
worked.  We  feel  it  is  reasonable  to  mention  two  of 
our  unsuccessful  attempts  so  that  these  practices  can 
be  avoided  in  the  future. 

Our  first  idea  was  to  use  the  water  as  the  support 
system  for  the  ice  plate.  However,  we  found  that  the 
forces  required  to  deform  the  plate  were  larger  than 
anticipated,  so  that  the  ice  plate  would  flood  when 
the  load  was  applied.  A  second  disadvantage  of  test¬ 
ing  the  ice  plate  while  in  contact  with  the  water  was 
that  a  strong  temperature  gradient  existed  across  the 
plate.  In  our  case  with  the  eoldroom  at  -10°C  and 
the  water  at  the  base  of  the  plate  at  +4°C,  a  tempera¬ 
ture  gradient  of  I4°C  existed  across  the  I  -  to  2-cm- 
thick  ice  plate.  This  temperature  gradient  compli¬ 
cated  the  characterization  of  the  experiment. 

The  second  method  that  we  used  with  limited  suc¬ 
cess  allowed  us  to  place  the  ice  plate  on  a  platform 
made  up  of  four  iron  bars  arranged  in  a  square  within 
the  tank.  These  bars  were  placed  into  the  tank  5  to 
10  cm  below  the  surface  of  the  water.  When  freezing 
of  the  ice  sheet  was  completed  and  the  ice  plate  had 
been  trimmed  to  size,  water  was  siphoned  from  the 
tank  until  the  water  surface  was  some  5  to  10  cm  be¬ 


low  the  top  of  the  support  bar.  The  ice  plate  came  to 
lest  on  the  support  structure  when  the  water  was  low¬ 
ered  below  the  level  of  the  constraining  bars.  The 
chief  drawback  with  this  method  was  that,  since  the 
bottom  of  the  plate  and  also  the  support  bars  were  wet. 
the  plate  tended  to  freeze  to  them.  This  complicated 
the  boundary  condition  of  the  plate  since  it  produced 
a  bending  resistance  at  the  plate  edge. 

The  method  that  finally  evolved  and  which  proved 
very  satisfactory  was  a  variation  of  the  method  de¬ 
scribed  above.  After  the  plate  was  cut  to  size  it  was 
removed  from  the  tank.  The  water  was  then  lowered 
and  the  plate  was  placed  on  the  supporting  structure 
for  testing  after  the  support  bars  and  underside  of  the 
ice  were  dry.  This  kept  the  plate  from  freezing  to  the 
bars  and  bending  resistance  at  the  boundary  was 
eliminated. 

For  future  work  we  feel  that  the  support  structure 
could  be  totally  removed  from  the  freezing  tank.  In 
this  manner  the  number  of  tests  could  be  increased 
because  a  second  plate  could  be  growing  while  the  first 
plate  was  tested.  This  would  increase  the  number  of 
plates  that  could  be  tested  in  a  given  time  interval. 

An  illustration  of  the  support  structure  for  the  ice 
plate  with  a  plate  in  position  and  the  load  applied  is 
shown  in  Figure  1 .  The  bars  that  support  the  plate 
are  32  mm  in  diameter  and  0.85  m  long.  The  plate  is 
supported  on  its  perimeter  by  the  bars  on  a  square 
that  is  0.77  m  on  each  side.  A  plan  view  of  the  sup¬ 
port  structure  including  the  transducer  location  is 
shown  in  Figure  2. 
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Figure  2.  Flan  view  of  the  support  structure 
and  ice  plate  showing  the  position  of  the 
transducers. 


To  facilitate  the  recovery  of  the  transducers  and 
the  loading  weights  after  the  ice  failed,  a  plywood 
platform  was  placed  beneath  the  ice  plate  This  plat¬ 
form  allowed  recovery  of  the  fractured  pieces  of  ice 
for  examination  at  the  completion  of  the  test. 

Acoustic  emission  monitoring  system 

We  monitored  the  acoustic  emissions  from  the  ice 
as  it  deformed  with  an  Acoustic  Lmission  Technology 
Corporation  Model  204B  monitoring  system.  This 
system  is  typical  of  acoustic  emission  systems  used 
for  monitoring  stress  waves  generated  f  rom  various 
sources  in  materials.  For  these  tests  the  data  were 
processed  in  terms  of  the  total  number  of  acoustic 
events  detected  over  the  duration  of  the  test. 

We  assume  that  the  source  of  the  acoustic  emissions 
in  our  ice  is  the  formation  of  small  cracks  and  tnicro- 
cracks.  This  interpretation  is  consistent  with  the  cor¬ 
relation  between  acoustic  emissions  and  small  frac¬ 


tures  in  ice  observed  by  (.old  ( I 960).  In  our  experi¬ 
ments  the  threshold  at  which  acoustic  emissions  are 
detected  is  well  below  the  limit  of  visible  fractures 
in  ice. 

Figure  3  is  a  plot  of  the  calibration  curve  t'or  the 
transducer  as  a  function  of  frequency  in  Mil/.  I  or  the 
amplifiers  used  in  these  experiments  the  bandpass  be¬ 
tween  -3-dB  points  was  100  to  300  kHz.  The  trans¬ 
ducer  calibration  curve  (Fig.  3)  suggests  that  we  would 
detect  the  acoustic  emissions  at  the  resonance  peak 
with  a  sensitivity  of-70  (IB  (re  I  V/pbar )  at  just  over 
150  kHz.  For  our  test  the  amplifier  was  set  to  a  gain 
of  84  dB  and  the  threshold  detection  level  set  between 
0.059  V  and  1 .43  V.  These  system  parameters  indicate 
that,  for  a  signal  to  be  detected,  the  signal  at  the  trans¬ 
ducer  had  to  have  a  minimum  intensity  of  0.001  3  Pa 
for  the  0.59-V  threshold  setting  and  an  intensity  of 
0.0285  Pa  at  the  detector  for  the  1 ,43-V  threshold 
setting.  (See  Appendix  B  for  details  on  these  determin¬ 
ations.) 


Figure  3.  Acoustic  emission  transducer 
calibration  curve. 


Frequency  ( MHz ) 


4 


Two  methods  were  used  to  couple  the  acoustic 
emission  transducers  to  the  ice  plate,  hot  the  first 
test  the  transducer  was  frozen  directly  to  the  ice 
plate  in  the  position  indicated  in  Figure  5.  For  suc¬ 
ceeding  tests  the  transducer  was  interfaced  with  the 
ice  surface  using  a  silicone  high  vacuum  grease. 

To  test  the  system  response,  a  0.5-mm  pencil  lead 
was  broken  on  the  plate  at  a  position  diagonally  op¬ 
posite  the  acoustic  emission  transducer.  This  demon¬ 
strated  that  a  signal  of  this  amplitude  could  be  de¬ 
tected  at  a  distance  of  approximately  0.75  m  from 
the  transducer.  For  our  experiment  we  assumed  that 
most  of  the  signals  detected  at  the  K4-dB  gain  level 
occurred  in  the  ice  plate. 

Displacement  transducers  and  data  recording 

A  number  of  different  displacement  transducers 
were  used  in  this  experiment,  as  indicated  in  Figme  2. 
The  load  was  applied  at  the  centei  of  the  ice  plate  and 
a  direct  current  deformation  transducer  (IK'DT)  was 
mounted  at  the  center  line  of  the  plate  12  cm  from 
the  central  point  of  the  plate.  A  lincai  variable  differ¬ 
ential  transducer  (  LVDT)  was  also  used  to  measure 
the  plate  deflection  at  a  point  1 1  cm  to  the  left  of  the 


plate  center  and  22  cm  below  the  IK'DT  as  viewed  m 
Figure  2.  As  indicated  ill  this  figme,  we  attempted  to 
use  two  lineal  potentiometers  hut  found  that  these  did 
not  have  the  scnsilivits  to  measure  the  deflections 
observed. 

The  iransduceis  and  the  acoustic  emission  preampli¬ 
fier  were  supported  oil  a  beam  that  extended  across  die 
lest  tank.  A  photograph  ol  the  lest  system  including 
the  loading  weights  is  shown  in  Figure  4.  In  this  photo¬ 
graph  we  have  labeled  each  test  piece. 

The  data  from  the  l)CI)l  and  die  acoustic  emission 
system  were  recorded  on  a  strip  cliait  lecoidei  Data 
from  the  I  VDT.  the  three  thermocouples,  die  DfDl. 
and  the  acoustic  emission  monitoring  system  were  re¬ 
corded  on  a  data  logger  The  interval  between  data 
acquisition  periods  varied  depending  on  which  phase 
a  particular  test  was  in. 

Mechanical  loading  system 

All  tests  were  conducted  in  a  constant  load  config¬ 
uration.  Weights  were  applied  at  the  center  ot  the  ice 
plate  as  indicated  in  Figures  I  and  2  To  eliminate 
noise  generation  between  the  load  and  the  ice  plate  a 
cloth  mat  was  used  to  isolate  the  load.  This  helped 


Load 


Figure  4.  Photograph  of  the  test  setup  with  each  piece  of  test  equipment  indicated. 


Acoustic  Events 


ensure  that  no  extraneous  acoustic  emissions  were 
generated. 

ANALYSIS 

Nevel  (1966)  Iras  analyzed  tire  time-dependent  be¬ 
havior  of  an  ice  plate  of  infinite  extent.  However, 
due  to  the  dimensions  of  our  experimental  ice  plate, 
it  was  necessary  for  us  to  seek  a  solution  in  terms  of 
finite  plate  theory.  It  was  also  necessary  ,  in  view  of 
our  experimental  arrangement,  to  consider  the  load 
to  be  distributed  over  a  finite  extent  of  the  ice  plate. 

In  Appendix  A.  the  elastic  solution  for  the  deflec¬ 
tion  of  a  finite  ice  plate  with  a  distributed  load  is  de¬ 
veloped.  The  solution  is  presented  in  terms  of  the  de¬ 
flection.  curvatures  and  bending  moments  developed 
in  the  ice  plate. 

In  order  to  examine  the  time-dependent  deforma¬ 
tion  of  the  ice  plate,  we  have  solved  for  the  viscoelas¬ 
tic  lesponse  of  the  plate  by  applying  the  correspondence 
principle  ol  linear  viscoelasticity  The  solutions  for  the 
viscoelastic  deformation,  curvatures  and  bending  mo¬ 
ments  arc  developed  in  Appendix  ('.  To  obtain  this 
solution  we  assume  that  the  ice  behaves  as  a  Maxwell 
fluid.  In  Appendix  C  we  also  develop  the  solution 
when  the  weight  of  the  plate  is  considered. 

EXPERIMENTAL  RESULTS 

In  this  experiment  six  ice  plates  were  tested.  Each 
of  the  ice  plates  was  giown  as  described  in  the  section 
on  experimental  procedure.  The  ice  plates  generally 
took  24  to  36  hours  to  grow  and  were  between  1  and 
2  cm  thick.  The  results  varied  considerably  from  test 
to  test  even  though  the  growth  history  ol  each  ice 
plate  was  similar. 

The  first  ice  plate  that  was  tested  had  an  average 
ice  thickness  of  I  .9  cm.  A  load  ol  1%  N  was  applied 
io  tire  center  of  the  plate,  and  the  plate  failed  approxi¬ 
mately  45  minutes  after  the  load  application.  Unfor¬ 
tunately.  few  data  points  were  recorded  (Fig.  5a ) so 
that  the  estimation  of  Maxwell's  constants  could  not 
be  reliably  made.  The  second  ice  plate  that  was  tested 
had  an  average  ice  thickness  of  1 .5  cm.  After  the  plate 
was  formed,  a  load  of  49  N  was  applied  to  the  center. 
The  plate  deformation  was  monitored  at  one  point, 
near  the  center  of  the  plate  (sec  Fig.  21- 

Figure  5b  shows  that  the  deformation  rate  of  the 
second  plate  achieved  a  constant  value  after  loading, 
and  continued  to  deform  at  that  rate  until  the  49-N 
load  w„s  removed  some  20  hours  after  its  application. 
The  maximum  deflection  of  the  ice  plate,  as  measured 
by  the  DCDT,  was  0.46  mm  in  the  20-hour  time 
interval. 


Tlie  acoustic  emission  count  from  the  initial  loading 
phase  of  this  test  shows  little  organization.  That  is.  we 
observed  no  regular  pattern  of  acoustic  activity  for  this 
loading.  However,  there  were  several  periods  of  in¬ 
creased  emissions  -  perhaps  indicative  of  local  yielding 
at  points  throughout  the  plate.  For  this  test,  the 
acoustic  emission  sensitivity  was  much  higher  than 
for  the  succeeding  tests.  It  was  set  to  detect  signals 
above  a  threshold  of  0.75  V  or  0.01 5  Pa  calculated  in 
accordance  with  the  procedure  described  in  Appendix 
B. 

When  the  load  was  removed,  the  ice  plate  was  al¬ 
lowed  to  relax  for  a  period  of  about  two  hours.  The 
initial  recovery  of  the  ice  plate  occurred  within  one 
hour  of  unloading.  Little  or  no  recovery  was  observed 
after  this  time. 

The  ice  plate  was  then  loaded  and  unloaded  through 
several  cycles  until  the  plate's  eventual  failure  (Fig.  5b). 
Unlike  the  initial  phase  of  this  test,  the  acoustic  emis¬ 
sion  count  during  the  loading-unloading  cycles  showed 
a  rather  regular  pattern,  typically  at  a  rate  of  100 
counts  per  minute.  This  may  indicate  that  local,  iso¬ 
lated  yielding  occurs  only  during  the  initial  loading  of 
an  unstressed  ice  plate. 

The  third  ice  plate  tested  was  subjected  to  a  differ¬ 
ent  load  history  than  the  first  two  plates.  In  this  test 
we  intended  to  apply  a  250-N  load  in  increments  of 
about  40  N.  The  total  load  was  to  be  applied  within 
one  minute.  In  this  test  the  ice  plate  (with  a  mean 
thickness  of  1.7  cm)  failed  when  the  load  reached 
157  N.  Acoustic  emissions  were  observed  throughout 
the  loading  of  the  plate  until  the  time  it  failed. 

In  the  next  test  (Fig.  5c).  the  plate  was  again  pre- 
loaded  with  a  49-N  load.  For  this  test  the  threshold 
intensity  at  which  acoustic  emissions  were  to  be  de¬ 
tected  was  increased  from  0.01  5  Pa  to  0.035  Pa,  which 
was  2.3  times  greater  than  in  the  previously  described 
experiments. 

During  the  application  of  the  49-N  load,  some  350 
emissions  were  detected  at  this  threshold  level.  The 
plate  had  a  mean  thickness  of  0.9  cm. 

The  specimen  was  then  unloaded,  allowed  to  recover, 
and  reloaded  first  to  49  N  and  then  to  1 54  N.  The  ice 
plate  failed  when  the  105-N  increment  was  applied  to 
bring  the  load  to  154  N.  Not  many  acoustic  emissions 
were  recorded  during  this  test,  due  perhaps  to  the  in¬ 
creased  threshold  level  setting. 

Tlie  last  two  tests  (Fig,  5d.  5e)  were  conducted  using 
nearly  equal  loads.  The  second  of  these  tests  (Fig.  5e) 
was  initially  loaded  to  93  N  and  allowed  to  deform  for 
a  22-hour  period.  The  load  was  then  increased  to  236 
N.  The  plate  continued  to  deform  for  1 1  hours  under 
this  load,  until  the  test  was  terminated.  The  ice  plate 
had  a  mean  thickness  of  1 .5  cm.  The  acoustic  emission 
data  from  this  test  are  unreliable  due  to  poor  coupling 
between  the  transducer  and  the  ice  plate. 
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Table  I .  Summary  of  tests  conducted. 
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Figure  6.  Elements  oj  the  Maxwell  model  with  the  parameters 
identified  as  specified  in  Table  I. 


The  reduction  of  the  load  deformation  data  when 
considering  the  initial  phase  of  loading  gave  the  values 
of  the  ice  plate  material  parameters  summarized  in 
Table  1 .  The  weight  of  the  ice  plate  cannot  be  ne¬ 
glected  in  comparison  to  the  applied  load  and  there¬ 
fore  has  been  taken  into  account  when  determining 
the  ice  plate  response  (see  Appendix  C).  The  Maxwell 
rheological  model  (Fig.  6)  material  parameters  ql/p{ 
(elastic  modulus)  and  p„  (viscous  time  constant)  were 
determined  and  are  reported  in  Table  1 .  Consider¬ 
able  scatter  in  the  obtained  values  is  probably  due  to 
the  grain  size  and  ice  plate  thickness  variation.  For 
instance,  the  ice  plate  of  test  3.  with  about  the  same 
thickness  as  test  I  and  a  smaller  applied  load,  failed 
instantaneously,  whereas  some  45  minutes  elapsed  be¬ 
fore  the  plate  of  test  1  collapsed.  We  found  that,  upon 
adding  a  143-N  increment  of  load  in  test  6,  the  result¬ 
ing  tangent  modulus  increased  from  9.04x  10s  to 
2.1 34 x  10s  N/m2  and  the  viscous  time  constant  de¬ 
creased  from  28.261  to  3.692  hours.  This  result  in¬ 
dicates  that  in  this  instance  subjecting  the  plate  to  a 
prior  load  (93  N)  produced  a  more  rigid  ice  plate. 

THIN  SECTION  ANALYSIS 

After  each  test  was  completed,  samples  of  the  ice 
were  collected  and  placed  in  a  cold  storage  area  at  a 
temperature  of  -30°C.  When  all  the  tests  were  com¬ 
pleted.  thin  sections  from  each  ice  plate  were  prepared. 


The  thin  section  analysis  indicated  that  there  was  a 
large  variation  in  the  grain  size  of  the  ice  crystals  in  the 
ice  plates.  (Seneially.  however,  the  ice  crystals  in  a 
single  ice  plate  were  fairly  consistent.  Single  crystals 
extended  the  full  depth  of  each  plate. 

Figure  7  shows  the  characteristic  grain  structure  of 
the  ice  plates  tested.  The  average  grain  diametei  is  de¬ 
termined  by  counting  the  numbei  of  grains  along  a 
length  of  straight  line  in  the  thin  section  (I  ig.  7b)  and 
dividing  that  number  into  the  length  of  line  over  which 
the  grains  are  counted.  The  last  column  of  Table  I 
lists  the  mean  grain  diameters  from  each  plate  found 
using  this  technique.  For  the  six  tests  conducted,  the 
mean  grain  diameter  varied  by  as  much  as  a  factor  of 
two  from  the  smallest  grains  to  the  largest  grains. 

The  stress  fringes  visible  in  Figure  7  arc  the  result 
of  thermal  stressing  introduced  during  thin  section 
preparation  and  are  not  associated  with  the  deforma¬ 
tion  of  the  ice  plate. 

SUMMARY  AND  DISCUSSION 

As  previously  indicated,  this  report  represents  an 
initial  investigation  into  the  flexural  properties  of  ice. 

It  was  also  intended  to  examine  the  feasibility  of  using 
acoustic  emission  techniques  for  these  investigations. 
Considering  the  limited  time  in  which  the  work  was 
carried  out,  we  feel  that  the  investigation  was  success¬ 
ful  in  meeting  these  objectives. 
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a  Surface  parallel  thin  seer  ion. 


h.  Surface  perpendicular  thin  section. 

Figure  7.  Thin  section  from  the  ice  plate  from  test  I  made  after  failure  of  the  ice  plate. 


One  of  (he  most  interesting  phenomena  we  ob¬ 
served  in  these  tests  had  to  do  with  the  failure  of  the 
ice  plates.  In  the  cases  in  which  the  ice  plate  failed 
at  some  time  after  application  of  the  load,  we  could 
not  detect  any  accelerating  creep  that  might  indicate 
an  impending  failure.  Similarly,  although  we  generally 


detected  acoustic  emissions  associated  with  loads  that 
produced  failure,  we  did  not  see  any  acceleration  of 
the  emission  rate.  This  finding  appears  to  indicate 
that  it  may  be  difficult  to  predict  an  impending  fail¬ 
ure  in  flexure.  However,  it  may  be  that  failure  will 
always  follow  when  the  deformation  rate  exceeds 


sunn-  value  or  when  acoustic  emissions  are  detected. 
These  points  should  be  investigated  in  lutuic  work. 
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APPENDIX  A:  ICE  PLATE  LINEAR  ELASTIC  RESPONSE 

The  linear  elastic  response  of  a  simply  supported  rectangular  plate  subjected  to  uniformly  dis¬ 
tributed  loading  (perpendicular  to  the  plane  of  the  plate)  over  part  of  the  plate  surface  (Fig.  A 1 )  is 
developed.  Small  displacement  theory  is  assumed.  The  Ray leigh-R itz  procedure  is  utilized. 

The  bending  strain  energy  U  stored  in  the  plate  during  deformation  (Sechicr  1952)  is 

{^ir 

where  r?  =  the  plate  deflection 
p  =  Poisson's  ratio 

D=  the  plate  flexural  rigidity  |/:7i,/12(l-pr  )|  (A2 

h  =  plate  thickness 
/■.'=  Young's  modulus. 

The  change  in  potential  energy  ll1  of  the  applied  load  per  unit  area  p  is 

a»c  h  +  d 
?  2 

-//  17  dxdy.  (Ac 


Since  simple  support  conditions  prevail  on  all  four  edges  of  the  plate,  the  deflection  is  expanded 
in  the  following  series: 


__  nmx  .  irnv 

*  =  ^«sinT  sinT~ 


Tlie  bending  strain  energy  and  change  in  potential  energy  of  the  applied  load  expression  are  reduced 
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Minimizing  the  total  energy  in  the  system  U-W  with  respect  to  the  amplitude  term  T?mn  yields: 

.  mite  .  nttd 
,,  sin  -r—  sin  -Tr¬ 
ibp  2a  2b 

where  m  and  n  are  odd.  The  middle  surface  curvature  changes  k  ,  k  ,  k  are  given  by 


r  igurc  A 1.  Rectangular  plate  subjected  to  distributed  loading. 


The  bending  moments  Mx,  A/y,  Mxy  are  given  by 
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liquation  A6  lor  the  deflection  amplitude  term  T) mn  and  eq  A4,  A7  and  A8  for  determining  the  de 
flection,  curvature  changes  and  bending  moments  constitute  the  elastic  solution  for  the  deflection 
of  a  finite  plate  subjected  to  a  distributed  load. 


APPENDIX  B:  ACOUSTIC  EMISSION  SYSTEM  SENSITIVITY 

The  sensitivity  of' acoustic  emission  systems  is  often  given  in  terms  of  the  gain  of  the  amplifica¬ 
tion  system  used  to  monitoi  the  acoustic  emissions.  This  gain  figure  (in  decibels)  is  generally  accep¬ 
ted  hy  acoustic  emission  workers  in  terms  of  an  implied  assumption  that  the  band  width  of  the 
acoustic  emission  (Ah  1  system  and  the  transducer  sensitivity  are  within  familiar  limits.  However,  to 
be  somewhat  more  specific  we  can  relate  the  Ah  sensitivity  to  either  the  signal  voltage  at  the  pre¬ 
amplifier  input  or  to  tire  stress  wave  intensity  at  the  transducers. 

In  the  text  we  state  that  we  used  a  system  gain  of  84  dB  and  our  threshold  level  was  set  between 
0.059  and  1 .43  V.  These  threshold  settings  indicate  that  for  a  signal  to  be  recognized  as  an  acoustic 
event,  its  amplitude,  when  amplified,  must  exceed  these  threshold  voltage  settings. 

The  numerical  gain  associated  with  84  dB  can  be  determined  from  the  formula 

(;=wJli/2°  (Bl) 

For  84  dB.  the  numerical  gain  is  (l  =  15,849  which  allows  us  to  calculate  the  signal  amplitude 
from  the  transducer  to  be  between 

0.059  V 

=  3.723*10^ 


9.033*  I0’5  V 


or  approximately  between  4  and  90  pV.  According  to  the  manufacturer’s  calibration  data,  the  peak 
sensitivity  of  the  transducer  used  is  -70  dB  (re  1  V/p bar)  at  a  frequency  slightly  larger  than  1 50  kHz 
(see  Fig.  3).  If  we  assume  that  our  signal  amplitude  corresponds  to  the  peak  transducer  sensitivity, 
we  can  calculate  a  nominal  amplitude  of  the  stress  wave  at  the  transducer. 

Using  the  formula  for  amplitude  gain 


*0  =  *'r  10 


where  E0  =  output  voltage  per  unit  stress 

/•'r  =  reference  value  (i.e.  I  V/pbar  or  10  V/ Pa) 
dB  =  transducer  sensitivity  (i.e.  -70  dB). 

Then  substituting  the  gtven  values  in  cq  B2,  we  find  the  calibration  value  h.0  to  be 


os t-70)  =  3  1  f,2x  I0'4 


or  equivalently 


F/0  =  3.162*10-- 


since  1  pbar  =  0.1  Pa. 
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So  it'  the  threshold  sensitivity  is  set  such  that  the  transducer  output  voltage  of  4  pV  or  90  pV  is 
equired  to  register  an  event,  lire  corresponding  stress  intensity  thresholds  are  between 


4x  IQ'6  V 
3J  62x  10'3  V/Pa 


I.265x  I  O'3  Pa 


and 


90  x  I0‘fi  V 
3.1 62x  I0‘3  V/Pa 


28.46x  1 0‘3  Pa 


which  are  the  values  reported  in  the  text. 

Although  all  acoustic  emission  equipment  manufacturers  now  calibrate  their  transducers  by  refer¬ 
encing  them  to  the  I  V/p bar  setting  there  are  no  data  to  suggest  tlia*  this  value  is  either  absolute  or 
correct.  Indeed  one  manufacturer  (Dunegan  / Endevco)  suggests  that  the  current  calibration  method 
which  relates  output  voltage  to  stress  is  not  valid  at  al\  Dunegan /Lndevco  suggests  that  the  proper 
calibration  should  be  in  terms  of  velocity  [i.e.  dB  re  I  V/(m/s)) .  The  stress  figures  are  calculated 
above  so  that  we  can  have  a  basis  of  comparison. 
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APPENDIX  C  ICE  PLATE  LINEAR  VISCOELASTIC  RESPONSE 


The  assumption  is  made  that  tee  is  incompressible  under  hydrostatic  stress  and  obeys  Maxwell’: 
model  (Fig.  h|  tor  devjaioric  stress  and  strain  (Nevel  1966).  The  viscoelastic  response  can  be  ob¬ 
tained  trom  the  corresponding  elastic  one  using  the  correspondence  principle  (Flugge  1975).  Thu 
making  the  following  substitutions 


p 

into  equation  Ah.  the  Laplace  transform  of  the  amplitude  term  r;mn  is 

, 1  +/V 

7  mu  “  ^  •» 


(Cla) 

(C!b) 

(C2a) 


where  ni  and  n  aie  odd  and 


nine  mrj 
sin  — _  silt  — — 
9h/i  2a  lb 


nnJ—2  +  —  V 

V2 


and  a  step  functitm  loading 
p(t)  =  pMt) 


(C2b) 


(C  2c) 


has  been  assumed.  Taking  the  inverse  Laplace  transform,  the  deflection  amplitude  term  r?mn  is 
given  by 


»Imn  =  Clt  +  p^il, 


(C3) 


Equation  A4  together  with  eq  C3  provides  the  plate  deflection  for  the  assumed  viscoelastic  material. 

The  middle  surface  curvature  changes  k%,  ky<  kxy  can  be  obtained  from  eq  A7  provided  that 
eq  C3  is  used  for  Tjmn .  Finally  the  bending  moments  Mx,  My,  Mxy  are  given  by 


M 


.  nny 
sin—** 
o 


n  i  nt  1  » /  i  ttlTlX  M1V 

—  +  r-  —  I  M 1 1  sin -  sin  —f- 

'  a  b 


XT'  T-'  2/|3  rmn  A/.t  nnx  fEZ 

w«y  =  ZL-* n*|C  *MneM-r  cos_r 


(C4  a) 


(C4b) 


(C4c) 


Solution  considering  the  weight  of  the  ice  plate 

f  or  simplicity  assume  a  -  b  and  c  =  d.  Then  eq  C2b  reduces  to 


.  mire  me 

0 ,  4  sin  — —  sin  — — 

,  _  96  pa_  2 a  2a 


7r6  a.h3 


timim2  +  n2  )2 


Deflections  r?mn  due  to  the  load  of  magnitude/?  extended  over  the  area  c  hv d  are  given  by  eq  03 
with  C  from  eq  05.  Deflections  JJmn  due  to  ice  plate  weiglit  of  magnitude  pt  extended  over  the  area 
a  by  b  are  given  by  eq  03  with  from  eq  C6: 


_  96  Pla  I 

C  ny  "  '  "" 

mi  »»ui  n6  q,h3  nm(m2  +n2  )2 


Keeping  only  the  first  term  in  the  series  »»  =  n  =  I .  front  eq  C5 


,,  24  pa  .  ?  TC 

C =  —  - -  sin2  — 

rr6  q./j3 


from  er|  (  6 


24  fV 
rr6  q.  /,3 


to  take  into  account  the  weiglit  of  the  ice  plate  multiply  C  by  1  +  C|w/C. 

Solution  considering  an  increment  in  applied  pressure 
Ap  at  time  t&p 

liquation  C2c  now  becomes 


p(t)  =  pA(/)  +  (Ap)  A(/WAp) 
and  its  Laplace  transform  is 


thus  eq  C3a  now  becomes 


.  .  '.-ft, 


C(1  +  ^V'^') 


(1  +/’,S) 


and  eq  C3  now  becomes 


nmn  -  ('[t+P,  Mt)\+C^  |(/-/ip)+p1A(/-/ip)|. 


Solufion  considering  an  increment  in  applied  pressure  Ap 
at  time  tAp  plus  the  weight  of  the  ice  plate 

Equation  C2c  now  becomes 

PU)=  (/>  +  />,  >.A<f  >+<-V>)A(/-/^p); 
ilms  eq  C  da  now  becomes  for  »>  =  ti  =  | 


ifc,=  ripl/l+e"^’^^  -!_ 

l  "  "  sin^- 

\  2a  i 


where 


rin.  .  24  />a4  ,  7T c 

C(P)  -  - - sm'  — 

2a 


qh ' 


l.quution  C3  now  becomes 
ri, ,  =  C 


♦  *  _L_ 

”  sin^ 
2a  J 


(C9) 


(CIO) 


(Cll) 


(Cl  2a) 


(Cl  2b) 


(Cl  3) 
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APPENDIX  D:  EQUIPMENT  LIST 


Instrument 

Mcimtlai  Hirer 

Model 

Quantity 

General  equipment  list 

Datalogger 

Monitor  l  abs. 

4. tOO 

1 

Regulated  power  supply 

Strip  chart  recorder 

Power  Mate  t  o. 

Soltee.  Inc. 

VP-67225 

l  inear  position  differential 
transducer  ( potentiometer) 

Bourns  Instr. 

l>\ 2001041  701 

T 

1  mear  voltage  ditterential 
transformer  1 1  \  1)1  1 

Schaevit/  lug. 

1  R-I00R 

LVDl  transducer  readout 

Schaevit/  1  ng. 

IR  100R 

Direct  current  differential 

transducer  (  D('l)  1  1 

Sanbourn  Man. 

7|X  DT-2000 

1 

Deadbolt  engine  heaters 

1  reeman 

S 

Variable  transformer 

Sup.  lice.  Co. 

1  lb 

1 

0  OOl-nim  calibration 

instrument 

Starrct  Co 

1 

1  hcrmocouples 

Salinometer  and  pipettes* 

Beckman 

RB2.(TL-G20 

4 

1 

Acoustic  emission  (AE)  equipment  list 
Battery  charger 

A  1  .  1  cell.  Corp. 

204  Bat.  Chg. 

1 

Preamp 

AT.  Tech.  Corp. 

0204010 

1 

Amplifier 

AT.  Tech.  Corp. 

0204022 

1 

AT  transducer 

AT.  Tech.  Corp. 

AC  1751. 

1 

*  I  or  saltwater  ice  plate  experiments. 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 


Xirouchakis,  P.C. 

Investigation  of  the  acoustic  emission  and  deforma¬ 
tion  response  of  finite  ice  plates  /  bv  P.C.  Xirouc.iak 
M.  Chaplin  and  \'.F.  St.  Lawrence.  danover,  N.tl.:  V.S. 
Cold  Regions  Research  and  Engineering  Laboratory; 
Springfield,  Va.  ;  available  from  .National  Technical 
Information  Service,  iy81. 

iv,  2i  p.,  illas.  ;  2d  cm.  (  CRREL  Report  81-6.  ) 
Prepared  for  Directorate  of  Military  Programs, 

Office  of  the  Chief  of  Engineers  by  Corps  of  Engineers 
U.S.  Armv  Coll  Regions  Research  and  Engineering  Lab¬ 
oratory  under  DA  Project  4A161132AT24 . 
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